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Abstract 

We show that a good fitting to the first two years of COBE-DMR observations 
of the two-point angular correlation function of CBR temperature is given by mod- 
els with a non-zero infrared cutoff kmin in the spectrum of the primordial density 
perturbations. If this cutoff comes from the finiteness of the universe, say, a topo- 
logical T3 model, we find kmin ~ (0.3 — l.l)7rHo/c with confidence level 95%. Such 
a non-zero kmin universe would also give a better match to the observations both 
of the RMS quadrupole anisotropy of CBR and of galaxy clustering. 
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Almost half a century ago, Infeld and Scliild ^ had pointed out that if the size of the 
universe is finite, the infrared cutoff k^nm is then non-zero, and the infrared divergences 
in quantum field theory are automatically excluded. However, it was shown later that 
the problem of infrared divergence in QED can be solved even in an infinite universe 
{kmin = 0), because the infrared divergence in the radiative correction can be precisely 
eliminated by adding the soft photon bremsstrahlung to the vertex correction 0. In this 
treatment we have, in fact, assumed kmin < AE/ch, where AE is the energy resolution 
of the detection system. 

For experiments done in physics laboratory the condition of kmin < AE/ch is always 
held. For instance, if the size of the universe is of the order of today's horizon cHq ^ ~ 
3000 h^^ Mpc (where h is the Hubble constant in unit of 100 km s~^Mpc~^), in order to 
detect soft photons with energy ~ kmin the energy resolution AE should be as small as 
HHq. The uncertainty principle requires then that the time needed to measure the cosmic 
infrared cutoff is ~ Hg ^. Therefore, it is impossible to detect the cosmic infrared cutoff 
kmin by local experiments [Q]. 

However, astrophysics do provide experiments, such as cosmic background radiation 
(CBR), which have lasted as long as the age of the universe, ~ Hq ^, and be able to measure 
kmin comparable with Hq/c. Using this idea, we proposed that the CBR anisotropy is 
an effective tool to detect the size of a small universe. We showed that the ratios of 
quadrupole moment of CBR anisotropy to higher order multipoles sensitively depend on 
the ratio L/(cHg^), where L is the size of the universe [Q. To compare this L/(cHo^)- 
dependence with the first year data of the COBE-DMR. [Q, it has been found that the 
lower limit to the size of a cubic T3 universe should be much larger than that given by 
the distributions of galaxies and QSOs 0. 

In this letter, we will show that a non-zero kmin is possible. We are motivated by the 
Two Years of COBE-DMIi observations of the CBR anisotropy (hereafter, the two-year 
data), which has been available recently 0. Compared with the first year data the 
quality of the two-year data on many aspects has been significantly improved. A new 
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result is found to be that the RMS quadrupole amphtude Qrms = 6 ± 3fiK is significantly 
less than the most likely quadrupole- normalized amplitude Qrms-ps^ which is 12.4l'^|;3/iK 
for n=1.6, or 17.4 ± 1.5^K for n = 1.0, where n is the index of the power-law spectrum 
of density perturbation. This Qrms-Qrms-PS difference can not be totally explained by 
cosmic variance ||^. On the other hand, smaller Qrms may indicate the lack of density 
fluctuations on the largest scales Therefore, it is valuable to study the possibility of 
explaining the Qrms-Qrms-PS difference by models with non-zero infrared cutoff kmin- 

Since the two-year data were reduced in the scheme of standard inflationary universe 
(f2 = 1, A = 0, and simply connected spatial hyper surf ace), we consider a Gaussian 
and adiabatic primordial perturbation with power-law spectrum in a flat universe. The 
large-scale fluctuations in the GBR temperature are given by : 

where y = {2cHq^, fl) is a vector of length 2cHq^ pointing to the direction d, on the sky. 
5(k) is the Fourier amplitude of the density contrast 5(r). For a power- law perturbation, 
(|5(k)p) can be written as 



(I wH = ^Ag^fc^ (2) 

where is a large rectangular volume, A; = |k| and Aq is a constant determined by the 
variance of the perturbed potential. 

The observed temperature fluctuations of GBR on the celestial sphere are usually ex- 
pressed by spherical harmonics AT/T{fl) = J2im(^T'^r^i^)j where Yi^{n) are the spher- 
ical harmonic functions. Defining a rotationally invariant coefficient Ci = J2m{WT\'^)j o'^s 
found from eqs.(l) and (2) that 

Q = ^'jy^ ^f-f E k-'jf{ky) , (3) 
where ji{x) is the spherical Bessel function. Since our defined Ci is dimensionless, the 
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COi?£^ quadrupole amplitude Q is related to our quadrupole amplitude C2 by Q = C2 T, 
where T is the mean temperature of GBR. 



If the universe is of kmin = 0, eq.(3) becomes ITT 



C^= H,iJ + ^) All k-~^jf{ky)dk, (4) 
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where the lower hmit of the integration in eq.(4) is taken to be zero , bcCBjUSG /^tti^tt, — . 
When n < 3, eq.(4) gives 

^0 _ 2/ + i ^o r(V)r(/ + ^) 

where is the quadrupole moment in a kmin = universe. 

Let us consider a kmin 7^ universe, say, a cubic T3 universe [D, which is constructed 
from a flat and infinite universe by the following identification on the 3-dimensional flat 
hypersurface: (xi + IL,X2 + mL,X3 + nL) = {xi,X2,X3) for all integers l,m,n. L is the 
size of the universe. In this case, the coefficients Ci should be directly calculated from 
eq.(3), i.e. 

Q = ^fJ\^^ Aly'-{ykm.nr E(%)""'jf (%) ' (6) 

where knim = 27r/L, and the summation covers all possible states of the wave vector: 
k = kmin{lif^i'^)- Obviously, Ci now depends on three parameters: 1) the amplitude Aq 
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of the perturbation (or the quadrupole moment C2 ), 2) the infrared cutoff k^in (or the 
size of the universe L), and 3) the index n of the perturbation spectrum. When fc^m 0, 
or L — > 00, one has Ci . Therefore, this fcmm- dependence of Ci provides an effective 
method to detect non-zero infrared cutoff kmin- 

In the COBE-T)MK observations, two measurements are independent: a) the two- 
point angular correlation function C{6) of CBR temperature, and b) the RMS quadrupole 
amplitude Qrms- The two-point angular correlation function C{6) is related to Ci by 

c{e) = Y,CiW\i)Pi{co^e), (7) 

where Pi{x) is the Legendre function, and W{1) is a window function, which is W{1) = 
exp{— 1/2[/(/ + /)/17.8^]} 1^. Comparing the measured correlation function to the model 
of eqs.(5) and (7), one can find a most likely quadrupole-normalized amplitude Qrms-ps- 



The result based upon the first year data does show Qrms-ps ~ Qrms 0- Therefore, the 
model of kmm = 0, i.e. eqs.(4) and (5), is consistent with the first year data. However, 
for the two year data, the difference between the RMS quadrupole and the most hkely 
normahzed quadrupole is as large as 2 — 4cr (Even taking the cosmic variance into account, 
the difference is still significant at 90% confidence level for = 1 Therefore, one would 
no longer be able to confidently say that the standard model (n = 1) is totally compatible 
with the current COBE-DMR observations. 

We use the standard technique to test the T3 models and to estimate the most likely 
model parameters (quadrupole amplitude C2 and size L) by a x^-minimization over the 
two-point angular correlation function. For a given n and y/L, we estimate C2 rms-ps by 
minimizing over the data: 

2 ^ [Cj - c{ei)f 

where Cj and cTj are, respectively, the observed values and errors of the angular correlation 
at 6i, and crcv{0) is the la cosmic variance of the C{6) |]T^. In eq.(8) we assumed that 



the variances in different bins are mutually independent. This approximation is probably 
suitable for our purpose, because it has been known, at least for the first year date, that 
the best fit amplitude of quadrupole does not significantly depend on the nondiagonal 
part of the covariance matrix of Ci |jl3|. Since acv is also proportional to C2 rms-PS, we 



adopt an iteration procedure to calculate C2 rms-PS- First we assume a zero and find 
out the best-fitting value of C2 rms-PS- Using this value we calculate the acv based on 
100 Monte Carlo realizations of Ci, and do the minimization again and find out a new 
fitting value of C2 rms-PS- Using this new C2 rms-PS we repeat the minimization and find 
another more accurate value of C2 rms-PS- The iteration procedure is stopped until the 
differences of and of C2 rms-PS between the two consecutive minimizations are less 
than 0.1%. The final C2 rms-PS and are our desired values. In fact the calculation 
converges very fast, and we need less than 5 iterations. Fig.l shows the goodness of 
the x^"fit of each model, i.e. the probability P{> Xmin) ^^at the experimental data are 
drawn from a realization of the model. It can be seen from Fig.l that for kmm ~ 



models, i.e. L y = 2cHq^, the P{> Xmin) of = 1.6 is much greater than that of 
n = 0.6 or 1.0. This result is the same as that of likelihood analysis done by the COBE 
group. Although the n = 1.6 case can be comfortably accommodated by the two-year 
data [acceptance probability P{> Xmin) ~ 40%], the smaller indices n < 1, which are 



favored more by current observational data of galaxy clustering [|I^], have much lower 
acceptance probability. 

Differing from the first year data, a remarkable feature shown in Fig.l is the high peak 
in P(> Xmin) space at y/L ~ 0.8. This means that the best fit value of the infrared cutoff 
is kmin ~ 0.8ttHq/c. In Fig. 2 we plot such a best-fitting curve as well as the two-year data 
of C{9). We also plot the best-fitting curve of an k^in = universe (n = 1.6) in Fig.2. 
It can be clearly seen from both Figs.l and 2 that considering the possible existence of a 
non-zero cosmic infrared-cutoff substantially improves the model fit to the experimental 
data (even for n = 1.6), and the most probable value of y/L is ~0.8, which is almost 
independent of n. At the 95% confidence level, we have 0.3 < y/L < 1.1 for n = 1, 
0.5 < y/L < 1.1 for n = 0.6 and y/L < 1 for n = 1.6. 
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Fig. 3 plots our best-fitting C2 rms-ps ^s a function of y/L, where the index n is taken 
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to be 0.6, 1.0 and 1.6 as well. We plot the measured value of C2 rms by the bold line 
in Fig. 3, and its la range by the dotted area. One can find that the models with L < 
l.2y = 2AcHq^ give good agreement between Qrms-ps and Qrms, i-e. their difference is 
no longer larger than 1 a. This result is also almost independent of the power law index 
n. 

Because a finite universe lacks the fluctuation power on scales larger than its size, we 
require a larger amplitude Aq of the density perturbation to fit with the observational C{6) 



|15|. However, larger Aq will lead to a stronger clustering on smaller scales. Therefore 
we should study whether the amplitude Aq in a kmin 7^ universe is compatible with 
observed clustering of galaxies on scales, say, 8 h~^Mpc. In order to calculate the density 
fluctuations on smaller scales from a given Aq perturbation, we need to choose the transfer 
function T{k) of linear growth The T{k) in turn is mainly determined by the matter 



6 



composition of the universe. Here we use the T{k) of the standard cold dark matter (CDM) 
model and of the cold plus hot dark matter (CHDM) model {^cdm = 0.7, = 0.3, and 
h = 0.5 [|16]). Fig. 4 presents the predicted values of ag, the rms density fluctuation of 
a sphere r = 8/i^^Mpc, as a function of y/L for n = 1. In Fig. 4 we also plot a result 
of o"8 given by statistics of galaxy distribution. By examining clusters of galaxies. White 
et al. suggest that in an Einstein-de Sitter universe is between 0.52 and 0.62. 
These limits are shown as the solid lines in Fig.4. Obviously, small universes of y/L ^ 1 
are not acceptable, because they produce excessive clustering on 8 h'^Mpc scale for both 
T{k). The models of y/L < 1 are acceptable when the CHDM T{k) is used. [Considering 
uncertainties in the estimated limits of erg and in the choice of T{k), we caution readers 
against placing overaccurate constraints of the as-test on the universe sizes]. 

With reasonable transfer functions, the best-fitting value of power index n = 1.6 for 
large universe (L ^ 2cHq^) is hard to reconcile with the observational data of galaxy 
clustering fl^. While the COBE-BMR data allows the finite universe of L ~ 1.2cHq^ 
to have n = 1, which can be brought into agreement with galaxy clustering for some 
reasonable transfer functions. Therefore, compared with models of kmin = (infinite) and 
large kmin (small) universes, the moderate-size (or kmin ~ ttHq/c) universe is in better 
agreement with the observations when both the COBE-DMR result and the clustering 
of galaxies are considered. 

In summary, in terms of the three independent and basic tests, i.e. C{9), Qrms and 
the galaxy clustering, a good survived model among those discussed in the paper is the 
one of a non-zero kmin universe. If this cutoff comes from the finiteness of the universe, 
such as a topological T3 model, the most probable value of kmin is 0.8ttHo/c. However, it 
does not mean that the kmin 7^ can only be given by a multiply connected topology like 
T3. There are other mechanism for kmin 0. For instance, in the standard inflationary 
universe, if the parameter Ntot in the inflation factor exp{Ntot) is ~ 54 (somewhat of fine- 
tuning parameter may be needed), one would also have an initial density perturbation 
with spectrum cutted off at kmin ~ ttHq/c. Regardless of these variousness of possible 
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explanations, one can conclude, at least, that the COBE-DMR. observation of AT/T is a 
powerful tool to probe the cosmic infrared cutoff. The first two years of the COBE-DMK 
observations reveal the possible existence of non-zero infrared cutoff in the spectrum of 
the primordial density fluctuations. It seems to be the first time to derive an "observed" 
values of the cosmic kmin with confidence higher than 95%. However, we should remember 
the difficulty in the measurement of Qrms due to foreground contaminations. Therefore, 
in order to discriminate various explanations, it is necessary to have other tests. The 



amplitude of octavopole C3 of AT/T would be one of such tests |T^. We beheve, more 
convinced results related to the non-zero cosmic infrared cutoff will be obtained as more 
observations of CBR AT/T become available. 
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Figure captions 

Fig.l P(> Xmin) ^ function of y/L. Here Xmm calculated by fitting the finite L 
universe model to the COBE-DMR two-point angular correlation function of the cosmic 
temperature fluctuations. The index n of power law perturbations is taken to be 0.6, 1.0 
and 1.6. The dot-dash line denotes P(> xLin) = 5%. 

Fig. 2 The observed angular correlation Ci{9)T^ and the best-fit curves by a) a T3 
model with n — 1 and y/L — 0.80 (solid line), and b) an infinite and flat universe with 
n — 1.6 (dashed line). 
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Fig. 3 The best fitted quadrupole, i.e. rms-P5) ^ function of the size of the universe. 
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The thick line denotes the RMS quadrupole measurement €2^^^, and the dotted area is 
its 1(7 region. 

Fig. 4 The predicted in a cubic T3 models. Two types of transfer functions are 
assumed: the standard CDM (dotted line) and the CHDM (dashed line). The solid lines 
are the upper and lower limits to erg, given by statistics of galaxy distribution [17]. 
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This figure "figl-l.png" is available in "png" format from: 
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